Spin-chirality decoupling in the one-dimensional Heisenberg spin glass with 

long-range power-law interactions 
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We study the issue of the spin-chirality decoupling/coupling in the ordering of the Heisenberg 
spin glass by performing large-scale Monte Carlo simulations on a one-dimensional Heisenberg spin- 
glass model with a long-range power-law interaction up to large system sizes. We find that the 
spin-chirality decoupling occurs for an intermediate range of the power-law exponent. Implications 
to the corresponding d-dimensional short-range model is discussed. 
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The issue of the spin-glass (SG) ordering has been 
studied quite extensively for years, and continued to give 
an impact on surrounding areas |l|-|3(. Nevertheless, the 
true nature of the ordering of SG magnets still remains 
elusive and controversial. Since the magnetic interac- 
tion in most of real SG materials is known to be nearly 
isotropic, they should be described as a first approxima- 
tion by the isotropic Heisenberg model. Recently, con- 
sensus appears among various numerical works that the 
isotropic Heisenberg SG in three dimensions (3D) ex- 
hibits a finite-temperature transition, while the nature 
of the transition still remains controversial 

It has been suggested in Ref.[4| that the model might 
exhibit an intriguing "spin-chirality decoupling" phe- 
nomenon, i.e., the chirality exhibits the glass order at a 
temperature higher than the standard SG order, Tqg > 
Tsg [5rLZ|- Chirality is a multispin variable representing 
the handedness of the noncollinear or noncoplanar spin 
structures induced by frustration. By contrast, Refs.[8|- 
Hfj| | claim that the 3D Heisenberg SG does not exhibit 
such a spin-chirality decoupling, only a single transition 
which is simultaneously SG and chiral-glass (CG). 

To get deeper insight into the behavior in physical di- 
mension d — 3, it is often useful to study the phenomena 
by extending the dimensionality to general d dimensions. 
In the limit of low d, the short-range (SR) Heisenberg 
SG exhibits only aT = transition in d = 1. In d = 2, 
recent calculations suggest that the vector SG model, ei- 
ther the three-component Heisenberg SG [ll[ or the two- 
component XY SG 12j, exhibits a T = transition ac- 
companied by the spin-chirality decoupling, i.e., the CG 
correlation-length exponent vqg is greater than the SG 
correlation- length exponent vsg- The spin-chirality de- 
coupling associated with a finite-temperature transition 
could occur, if any, in d > 3. In the opposite limit of 
high d, the SR Heisenberg SG model in infinite dimen- 
sions d — > oo reduces to the mean- field (MF) model, i.e., 
the Sherrington-Kirkpatrick (SK) model. The Heisen- 
berg SK model is known to exhibit only a single finite- 
temperature SG transition, with no spin-chirality decou- 
pling. In high but finite d, Monte Carlo (MC) result of 
Ref.[13| suggested that the spin-chirality decoupling did 
not occur in d = 5, but might occur in d — 4. Reflect- 
ing an intrinsic difficulty in thermalizing large systems in 



high dimensions, however, the true situation still remains 
largely unclear. 

In the present paper, we attack the issue of the spin- 
chirality coupling/decoupling in the Heisenberg SG from 
a different perspective. Namely, we study a different type 
of Heisenberg SG model, i.e., the one-dimensional (ID) 
Heisenberg SG with a long-range (LR) power-law inter- 
action proportional to l/r a (r is the spin distance). In 
the limit of sufficiently large a — > oo, the model reduces 
to the standard d = 1 model with a SR interaction. In 
the opposite limit of a — > 0, the model reduces to an 
infinite-range model, i.e., the SK model corresponding 
to d = oo. Hence, varying a of the ID LR model might 
be analogous to varying d in the SR model. Indeed, this 
correspondence was supported by recent studies by Katz- 
graber and Young |l4j and by Leuzzi et al [l5[ for the 
Ising SG. These authors have suggested more detailed 
correspondence between d of the SR model and a of the 
ID LR model, e.g., (i) the upper critical dimension d = 6 
corresponds to a = 2/3, (ii) the lower critical dimension, 
which lies between d = 2 and 3, corresponds to a = 1, 
and (iii) d = 3 corresponds to a ~ 0.9. 

Advantages of studying such ID models might be 
threefold. First, systems of large linear size L, never 
available in high dimensions, can be studied (up to 
L = 4096 in the present calculation). Second, one can 
continuously change and even fine-tune the parameter a 
playing the role of effective "dimensionality" , while it is 
impossible to continuously change the real dimensional- 
ity d in the SR model. Hence, by studying the properties 
of the ID model with varying a, one might get an over- 
all picture concerning how the "coupling vs. decoupling" 
behavior depends on the effective dimensionality. Third, 
certain analytical results based on the renormalization- 
group (RG) calculations are available in ID, which might 
be utilized in interpreting the numerical data. 

Indeed, RG calculations, though did not take account 
of the possibility of the spin-chirality decoupling, sug- 
gested that the model exhibited a rich ordering behav- 
ior with varying a 0, [I?}- For tr < 2/3, the Gaus- 
sian fixed point is stable and the model exhibited a 
finite-temperature SG transition of the MF type. For 
2/3<cr<l,a non-trivial LR fixed point becomes sta- 
ble leading to a finite-temperature SG transition char- 
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acterized by the non-MF exponents. In particular, the 
critical-point-decay exponent is determined solely by the 
power describing the spin- spin interaction, and is given 
by rj SG = 3 - 2cr [H, [lZj]. For a > 1, the SG transition 
occurs only at zero-temperature with rise = 1- 
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FIG. 2: The (inverse) size dependence of the crossing tem- 
peratures of £cg/L, £,sg/L, and xsg/L 2 ~ tisg for the case of 
a — 0.9 (a), and of a = 1.0 (b). Lines represent power-law 
fits with an exception of the SG data for a = 1.0 which is a 
logarithmic fit: See the text for details. 



FIG. 1: The correlation-length ratio versus the temperature 
for the chirality (a), and for the spin (b), for a — 0.9. The 
arrows indicate the bulk chiral-glass and spin-glass transition 
points. The inset of Fig.(b) represents the temperature de- 
pendence of the spin-glass susceptibility ratio, xsg/ L 2 ~ vsg 
with 2 - t]sg = 2a - 1 = 2 x 0.9 - 1 = 0.8. 

Meanwhile, it remains to be seen how the spin-chirality 
decoupling arises in this ID model with varying a. Since 
the MF Heisenberg SK model does not show the spin- 
chirality decoupling, the spin-chirality decoupling associ- 
ated with a finite-temperature transition should be real- 
ized, if any, only in the intermediate range of <r, near or 
below (7 = 1. Thus, we study here both the spin and the 
chiral orderings of the model by large-scale MC simula- 
tions, varying a in the range 0.7 < a < 1.1, which spans 
the non-MF regime. Our numerical results indicate that 
the model exhibits the spin-chirality decoupling in the 
range 0.8 IS o % 1.1, while the usual spin-chirality cou- 
pling behavior occurs for a ^ 0.8. 

The Hamiltonian is the ID classical Heisenberg model 
with a random LR power- law interaction Jy , 

H = — 22 JijSi ■ Sj , (1) 

<ij> 



where Si = (Sf, Sf , S? ) is a three-component unit vector 
at the i-th site, and the < ij > sum is taken over all spin 
pairs on the lattice once. The interaction decays with 
a geometric distance as a power-law, 

J v = C ^~' C = V ^ ( 2 ) 

ij 2^i<ij> r ij 

where is an independent random Gaussian variable 
with zero mean and standard deviation unity. Periodic 
boundary condition is applied by placing L spins on a 
ring. Then, the geometric distance between the spins at 
i and j is given by ry = {L/it) sin(7r \i — j\ /L). 

We perform extensive MC simulations for various val- 
ues of a in the range 0.7 < a < 1.1 (a p reliminary report 
for a = 1.1 was presented in Ref . [la] ) . We shall show 
below mainly the results for a = 0.9 and 1.0. The lat- 
tice sizes studied are L= 128, 256, 512, 1024, 2048, and 
also 4096 for some a. Sample average is take over 896 
(for L < 2048) and 256 (L = 4096) independent bond 
realizations for a = 0.9, while 896 (L < 1024) and 256 
(L = 2048) for a = 1.0. We use a single-spin- flip heat- 
bath and an over-relaxation method combined with the 
temperature-exchange technique. The over-relaxation 
sweeps are repeated 5 times per every heat-bath sweep, 
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which is set as our unit MC step. Equilibration is checked 
by monitoring: i) All the "replicas" travel back and 
forth many times (typically more than 10 times) along 
the temperature axis during the temperature-exchange 
process between maximum and minimum temperature 
points, whereas the relaxation due to single-spin flip is 
fast enough (both chiral and spin autocorrelation times 
about 20 MC steps or less) at the maximum tempera- 
ture: (if) All the measured quantities converge to stable 
values 0j- 

The local chirality at the i-th site Xi is defined for three 
neighboring spins by the scalar x% = Si+i ■ (Si x Si-i)- We 
simulate two independent copies of systems with identical 
interaction sets Jij subject to mutually different random- 
number sequences and spin initial conditions, and mea- 
sure fc-dependent overlaps both for the spin and for the 
chirality. The fc-dependent chiral overlap, q x (k), is de- 
fined as an overlap variable between the two replicas (1) 
and (2) by 
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From the chiral overlap, we calculate the CG suscepti- 
bility XCG via the second moment of its k — compo- 
nent, xcg — L[{\q x (Q)\ 2 )], where (•••) denotes a ther- 
mal average and [• • • ] an average over the bond disorder. 
Essentially the same definitions also apply to the spin 
except that an appropriate overlap becomes a tensor in 
spin space [7]. Finite-size correlation length of the ID 
LR model is then defined by 



1 ( [(g(Q) 2 )] 

2sin(fc m /2) U<9(*m) 2 >] 
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(4) 



for both SG and CG, where fc m = ^ 03- 

We begin with the case of a = 0.9 in the midst of 
the non-MF regime. In Fig.l, we show the correlation- 
length ratios for the chirality £,cg/L (a), and for the spin 
£,Sg/L (b). These quantities are dimensionless so that 
the data of different L should cross or merge at the ex- 
pected transition point asymptotically for large L. For 
smaller sizes, the crossing temperatures T cross (L) of the 
spin £sg/L and of the chiral £cg/L are almost com- 
mon, exhibiting only a weak L-dependence. For larger 
sizes, T cross (L) of the spin £,sg/L shift down to lower 
temperatures, while those of the chiral £cg/L remain 
to be nearly L- independent, deviating from T cross (L) of 
the spin £sg/L. In Fig. 2(a), we plot T cross (L) of the 
£cg / L curves between of the sizes L and 2L as a func- 
tion 1/L, together with the corresponding ones of the 
£,Sg/L curves. 

While the SG and CG susceptibilities are dimensionfull 
quantities, they can be made dimensionless by dividing 
them by L 2 ~ n where rj is a critical-point-decay exponent. 
Generally, the exponent r\ is not known in advance, but 
in the case of the present LR interaction, the SG expo- 
nent rjsG is given by rjsc = 3 — 2a. Making using of this 



property, we plot in the inset of Fig. 1(b) the temperature 
dependence of the SG susceptibility ratio xsg/L 2 ~ r i SG 
where 2 — ijsg = 2<t — 1 = 0.8. As can be seen from 
the figure, the data of different L exhibit a crossing be- 
havior as expected for the dimensionless quantity. For 
the CG susceptibility ratio, this type of analysis has only 
restricted utility because of the lack of our knowledge of 
the chiral rjcG value. 
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FIG. 3: The correlation-length ratio versus the temperature 
for the chirality (a), and for the spin (b), for a — 1.0. The 
arrow indicates the bulk chiral-glass transition point. The 
inset of Fig.(b) represents the temperature dependence of the 
spin-glass susceptibility ratio, \sg / l L 2 ~' nsG where 2 — t)sg = 
2cr-l = 2xl.0-l = 1.0. 

In order to estimate the bulk CG and SG transi- 
tion temperatures quantitatively, we need to extrapolate 
Tcross(L) to L = oo. Such an extrapolation is done here 
on the basis of the relation, T cross (L) — T cross (oo) cL~ e 
with 8 — v^ 1 + lu, where v and u are the correlation- 
length and the leading correction-to-scaling exponents, 
respectively, while c is a constant. For the SG, we per- 
form a combined fit of T cross (L) of both £sg/L and of 
Xsg/L 2 -' 1sg , where a common T cross (oo) and 9 are as- 
sumed. The standard x 2 -analysis then yields Tsg = 
0.086 ± 0.003 and 6 = 0.44 ± 0.07, with x 2 /DOF= 1.24 
and the associated fitting probability Q = 0.29. The 
smallness of the obtained error bar of T$g is due to 
the fact that the two independent T cross (L) are used 
in the fit, each approaching Tsg either from above or 
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FIG. 4: The a versus temperature phase diagram of the 
ID Heisenberg spin glass with a LR power-law interaction 
decaying with a distance r as oc r~ a . The red (blue) points 
are the chiral Tcg (the spin Tsg) transition temperature. The 
spin-chirality decoupling occurs in the range 0.8 5 cr < 1.1, 
while more standard coupling behavior occurs in a $ 0.8. 

below. For the CG, we have T cross (L) of £,cg/L on ly> 
which yields T CG = 0.105 ± 0.003 and = 1.2 ± 1.4 
(xVDOF^ 1 - 14 and Q = 0.32). The smallness of the er- 
ror bar of Tqg is due to the fact that T cross (L) of ^cg/L 
exhibits a nearly L-independent behavior. Hence, Tcg 
turns out to be higher than Tsg by about 20%, suggest- 
ing that the spin-chirality decoupling occurs for a = 0.9. 

On decreasing a from a = 0.9, one tends to approach 
the spin-chirality coupling regime. Indeed, for a = 0.8, 
we obtained via similar analyses (the data not shown 
here) T CG = 0.158 ± 0.008 and T SG = 0.159 ± 0.002, 
which suggests that the spin and the chirality might order 
simultaneously at a = 0.8. 

By contrast, on increasing a from a = 0.9, ones ap- 
proaches the Tsc — regime, with the upper-critical cr, 
a = 1. In Fig. 3, we show the correlation-length ratios 
£,cg/L (a) and £sg/L (b) for cr = 1.0, together with 
the SG susceptibility ratio xsgI L 2 ~ vsg . In contrast the 
cr = 0.9 case, Xsg/L 2 ~ vsg does not show crossing in 
the investigated T and L range. In Fig. 2(b), T cross (L) 
°f £,cg/L and of £,sg/L as well as those of xsG/L 2 ~ nsG 
are plotted as a function of 1/L. For the chirality, a 
fit of T cross {L) of icc/L yields T CG = 0.045™? and 



9 CG = 0.34±0.34 (x 2 /DOF=0.15 and Q = 0.70). Hence, 
Tec is likely to be nonzero at a — 1.0. For the spin, the 
decreasing tendency of T cross (L) with L becomes pro- 
nounced. In fact, a power- law fit becomes unstable here, 
leading to an indefinitely negative Tsrj-value. Rather, 
a logarithmic fit of the form T cross (L) — b(h\L + c)~ e , 
expected for the T = transition at the upper-critical 
a, yields an acceptable fit with 9 ~ 2.1 (x 2 /DOF= 3.88 
and Q = 0.0044) as shown in Fig. 2(b). This observation 
supports the T — SG transition theoretically expected. 

Similar analyses are made for other values of a includ- 
ing cr = 0.7,0.8,0.85,0.95,1.1. We then find the decou- 
pling behavior with Tcg > Tsg for a = 0.85,0.95, and 
the coupling behavior with Tcg — Tsc for cr = 0.7, 0.8 
(details will be reported elsewhere). The obtained re- 
sults are summarized in the a — T phase diagram in 
Fig. 4. We tend to have larger error bars toward the 
upper-critical cr = 1, because we have only one kind of 
Tcross(L) and the data exhibit significant downward cur- 
vature there. The spin-chirality decoupling occurs in the 
range 0.8 IS cr "S 1.1. By contrast, the standard spin- 
chirality coupling behavior Tsg = Tcg is realized for 
cr < 0.8. The para-CG phase boundary might go beyond 
cr = 1, touching the T-axis separately from the CG-SG 
phase boundary. 

Although the correspondence between a of the ID LR 
model and d of the SR model is by no means exact, it 
might be interesting to deduce on the basis of Fig. 4 the 
behavior of the c?-dimensional Heisenberg SG with SR 
interactions. As mentioned, recent theoretical analyses 
suggest that the physical dimension d — 3 corresponds 
to cr slightly smaller than unity. Hence, our present con- 
clusion that the spin-chirality decoupling occurs in the 
relevant range of cr, 0.8 5 cr < 1.1, gives indirect sup- 
port to the spin-chirality decoupling occurring in d = 3 
in the SR Heisenberg SG. Anyway, the phase diagram of 
Fig. 4 serves to grasp an overall ordering behavior of the 
Heisenberg SG from a wider perspective. 
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